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RESUMEN 

Analizamos la distribucion del material interestelar en la vecindad de las es- 
trellas de tipo O: HD 38666, HD 124979, HD 163758 y HD 171589, a fin de 
investigar la existencia de burbujas interestelares asociadas a las estrellas. La 
localizacion de las estrellas lejos del piano galactico, favorece la formacion de 
burbujas de gran extension. En base a datos del Hi, continuo de radio y CO, 
e imagenes en el infrarrojo lejano que revelan la emision del polvo interestelar, 
investigamos la distribucion del gas neutro e ionizado en estas regiones. Re- 
portamos en este trabajo la deteccion de cavidades y envolturas en expansion 
en hidrogeno neutro, asociadas a las cuatro estrellas masivas. Se discute la 
posibilidad de que estas estructuras hayan sido originadas por la accion de los 
vientos estelares sobre el gas circundante. 

ABSTRACT 

In our search for interstellar bubbles around massive stars we analyze the en- 
virons of the 0-type stars HD 38666, HD 124979, HD 163758, and HD 171589. 
The location of the stars, which are placed far from the galactic plane, favors 
the formation of large wind bubbles. We investigate the distribution of the 
neutral and ionized gas based on Hi, CO, and radio continuum data, and that 
of the interstellar dust based on far infrared IRIS images. Here we report 
the discovery of neutral gas cavities and slowly expanding shells associated 
with the four massive stars. IR and optical counterparts were also detected for 
some of the stars. We discuss the probability that the features have originated 
in the action of the stellar winds on the surrounding gas. 

Key Words: ISM: bubbles — stars: individual (HD 38666, HD 124979, HD 163758, 
HD 171589) 



1. INTRODUCTION 

The strong stellar winds of massive stars interact with their environs creat- 
ing interstellar bubbles that are detected in a large range of wavelengths. UV 
photons with energies \\v > 13.6 eV ionize the inner p art of the expa nding bub- 



the r adio continuum ranges as shell-like thermal sources (jCoss fc Lozinskava 



|xp; 

bles, which are generally detected both in the optical (iLozinskaval 1982) and in 



19951 ). Should the ionizing front be trapped within the expanding envelope. 



iJnstituto Argentino de Radioastronomi'a (CCT-La Plata, CONICET), C.C. 5, 1894, 
Villa Elisa, Bs. As., Argentina 
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interstellar bubbles have an outer neutral region tha t can be identified both 
in the Hi 21-cni line emission and in molecular lines ( Cappa et al. 20051 ) 

Many interstell ar bubbles are detected neither in the radio continuum nor 

in op tical lines (e.g. Cappa fc Benaghall998 : Arnal et al. ICichowolski fc Arnall 
20041 ). A low ambient density may be responsible for the undetected radio 



continuum or opt ical emission, as was suggested for bubbles in the LMC by 



Naze et al.l (|2002l ). As a consequence, the analysis of the Hi emission distri- 
bution is an important tool to investigate the presence of such structures and 
allows discrimination in distance based on galactic rotation models. 

Data at radio and infrared wavelengths offer an opportunity to investigate 
the characteristics and distribution of the ionized and neutral gas associated 
with these bubbles, as well as those of the interstellar dust. 

A low-velocity massive star with constant stellar wind parameters located 
in an homogeneous interstellar medium with uniform density creates a spher- 
ical stellar wind bubble. Aspherical bubbles appear as the result of relaxing 
some of these hypothesis. The combination of a powerful ste llar wind and 



a high space velocity makes this scenario changes drastically. IWeaver et al 



(|1977l ) pointed out that a massive star moving supersonically with respect to 
the surrounding gas originates an aspherical stellar bubble, elongated in the 
direction of the movement of the star. A w ind bow-shock can develop ahead 
of the star (|Wilkinlll99i iRaga et al.lll997h . As Ivan Buren fc McCravl |l988h 
have shown, these shocks can be identified in the far IR emission. 

Here we present a large scale study of the interstellar medium around four 
southern 0-type stars: HD 38666, HD 124979, HD 163758, and HD 171589, 
the first two of them classified as runaway stars. Our aim is to investigate 
the action of the UV photon flux and the stellar winds of these stars on their 
surroundings. The selected stars are located at |&| > 3°, allowing for the 
formation of relatively large structures which are easily identified using low 
angular resolution data. 

In the following sections, we analyze the interstellar medium in the envi- 
rons of the stars looking for cavities and shells in the neutral gas that have 
originated in the action of the massive stars on their surroundings. These 
studies are important to improve theoretical models on the interplay between 
massive stars and the surrounding gas. 



2. DATA BASES 



The Hi 21-cm line data analyzed in this paper were extracted f rom the Lei- 
den A rgentine Bonn Survey of Galactic Hi ( Kalberla et al. 20051 : Arnal et al. 



20001 ) observed using the radiotelescope of the Instituto Argentino de Radioas- 
tronomfa. The 21-cm data span the velocity range from -250 to -f250 km 
and were obtained with a velocity resolution of 1.3 km s~^. The angular 
resolution is 30'. 

Radio continuum da ta at 4.85 GHz from the Parkes-MIT-NRAO (PMN) 
Southern Radio Survey (ICondon et al.lll993l) are available on lv for HD 124979 
and HD 171589, while data at 408 MHz |Haslam et al.|[l98i ) and at 35 MHz 
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TABLE 1 

MAIN PARAMETERS OF THE 0-TYPE STARS 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


Name 


(l,b) 

n 


Spectral 
Classification 


d 

(kpc) 


V 
(mag) 


(B-V) 
(mag) 


(]" 

^spc 

(kpc) 


cosb^ 
(mas yr""*^) 


Mb 

(mas yr~^) 


z 

(pc) 



HD 38666 


237.29,-27.10 


09.5 V" 


0.53" 


5.17" 


-0.27" 


0.8 


+23.48±1.1 


-3.52±1.1 


-255 


HD 124979 


316.40,+9.Q8 


08 HI ((f))'' 




8.52" 


1.05" 


3.7 


-8.8±1.1 


+8.2±1.1 


+590 


HD 163758 


355.36,-6.10 


06.5 la (f)" 




7.32" 


0.03" 


3.5 


-2.5±1.5 


-5.2±1.5 


-374 


HD 171589 


18.65,-3.09 


07 H (f)" 


1.15^-1.5'' 


8.29" 


0.29" 


3.0 


+6.0±1.5 


-4.1±1.5 


+160 



"COS catalogue (|Mai'z-Apellaniz et all (|200J)) 
''See text 



IKozokl (| 19851) 

• iGarmanv et all (|l982h 

°Spectropliotometric distances derived in this paper. See text for details 
•^Derived from Tycho-2 catalogue (|Hog fc Murdinlboooh . 



( Dwarakanath &: Udava Shankai IQQOl ) are available for the four target stars. 
The angular resolutfons are 5', 0?85, and 26' x 42'/cos((5-14°), for the surveys 
at 4.85 GHz, 408 MHz, and 35 MHz, respectively. 

The analysis of the molecular gas distributi on for the reg i ons of HD 163758 
and HD 171589 was based on the CO survey bv lDame et al.l (|200l[ ). which has 
an angular resolution of 8'.8. The velocity coverage and the velocity resolution 
are -160 to +160 km s"'^ and 1.3 km s^^, respectively. The rms noise level is 
0.2 K. 

The dust distribution in the region was analyzed using the IRI^ database. 
These images are a new generation of IRAS images that benefit from a better 
zodiacal light subtraction, calibration, a nd a better destriping. IRIS image s 
have an angular resolution of around 4' ( Miville-Deschenes Sz Lagache 20051) . 

The optical im ages were obtaine d from the Full- Sky H- Alpha survey (H- 
Alpha Composite. iFinkbeineri (|2003l )) 



3. TARGET STARS 

Table 1 summarizes the stellar parameters relevant to this study: the 
name of the stars and their galactic coordinates are listed in the first two 
colu mns. The spectral types, in dicated in the third column, were obtained 
from iMaiz-Apellaniz et al.l (|2004( )(GOS catalogue). Distance estimates taken 
from the literature are given in column 4. Columns 5 and 6 list the visual 
magnitudes V and color indices {B — V), respectively. W e used these values 



together with absolute magnitudes from lVacca et al.l (|l996r ) and intrinsic color 



'http:/ /www. ias.fr/IRIS 
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indices from IWegnen (|l994l ) to derive the spectrophotometric distance listed 
in column 7. The components of the proper motion in the galactic coordinate 
system listed in colum ns 8 and 9 were calculated from the Tycho-2 catalogue 
(jHog fc Murdinll2000l ). Finally, z-distances to the galactic plane are shown in 
col. 10. They were derived by adopting the distances listed in col. 7. 

Based on its large proper motion, HD 38666 (/i Col) was identified as a 
runaway star belonging to OriOBl (|Blaauwi,1961i) . 



HD 124979 has been classified as an 8 ((f)) sta r by 



MacConnell k, Bidelman 



( 1976 ). Following th e Of typificat i on by Walborn I l97l[ ). we have adopted a 



luminosity class III. iMason et al.l (jl998t ) analyzed the binarity status, and 
concluded that it maybe a spectroscopic binary, tagging it as "SBl?". The 
measured proper motion and radial velocity (—68 km s~^) allow to classify 
this star as a runaway. 

No information about the interstellar medium in the environs of the target 
stars was found in the literature. 



4. Hi emission TOWARDS THE SELECTED STARS 

Figure 1 {a-d) displays the Hi 21-cm line profiles towards the 0-type stars. 
They were obtained by averaging the Hi emission within an area of about 
3°x3° centered at the position of the massive stars. Our aim is to analyze 
the main characteristics of the neutral gas emission towards the selected stars 
and its relation to the galactic spiral structure. 

The Hi emission profile in the region of HD 38666 (Fig. la) reveals a weak 
Hi gas component centered at v « km s"^ (all velocities in this paper are 
referred to the LSR). Because of the high galactic latitude of the star, this 
material corresponds to local gas. Bearing in mind the distances to HD 38666 
listed in Table 1, we expect that gas related to this star has low positive 
velocities. 

The Hi emission profile in the region of HD 124979 (Fig. lb) reveals three 
components centered at v « -2, v « -20, and v « -50 km . The for- 
mer gas component is associated with th e local region. The circular galactic 
rotation model by Brand &: Blit j (1993) predicts near and far kinematical 
distances c? « 1.4 and 11 kpc, respectively, for the component at -20 km 
s"""^, while for the component at -50 km s""'^, d « 3.7 and 8.5 kpc, respec- 
tively. The near kinematical distances are compatibl e with the lo cations of 
the Sagittarius-Carina an d Scutum-Crux spiral arms ( Russeill 2003 ). Accord- 
ing to the analytical fit bv lBrand fc Blita . gas in the environs of this star, i.e. 
at 3.7 kpc should have velocities of about -50 km s~^. 

In the line of sight to HD 163758 (Fig. Ic), Hi emission was detected from 
-40 to +40 km^^. The profile shows a main peak of neutral gas emission 
centered at v w -1-5 km s~^, probably related to the local spiral arm. The fact 
that positive velocit ies are forbidden in the fourth galactic quadrant within 
the solar circle (e.g. Brand fc Blitz 1993h indica tes that mos t of this material 
is local and associated with the Gould's belt (jOland 119821 ). Gas linked to 
this star at a distance of 3.5 kpc should have velocities of about -10 km . 
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Fig. 1. Hi profiles obtained by averaging the Hi emission in regions of 3° x3° around 
the O-type stars, (a) HD 38666; (b) HD 124979; (c) HD 163758; and (d) HD 171589. 
The intensity scale is given in brightness temperature. Velocities are referred to the 
LSR. 
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Note however that kincmatical distances have large uncertainties in this hne 
of sight. 

The Hi emission towards HD 171589 (Fig. Id) was detected within the 
velocity interval -40 to +100 km s~^. The neutral gas emission shows major 
peaks at v w 0, V +12 km s^^, and v « +40 km s~^. The lower velocity peak 
belongs to the local spiral arm, while peaks at positive velocities correspond 
to near and far kinematical distances d « 1.3 and 15 kpc, and d « 3.4 and 
12.5 kpc, respectively. These gas components are probably related to the 
Sagittarius-Carina and Scutum-Crux spiral arms, respectively (Russeil 2003). 
A relatively low emiss ion component can also be identified at « -25 km s~^. 
The analytical fit by iBrand &: Blita predicts that gas with this velocity is 
located at « 20 kpc, well beyond the solar circle. Material related to the star, 
placed at d < 3 kpc, should appear with velocities lower than +33 km s^^. 



5. Hi structures AND THEIR COUNTERPARTS AT OTHER 

FREQUENCIES 

In order to look for Hi structures linked to the selected stars, we analyzed 
the neutral hydrogen distribution in their environs. A series of (/, b) maps at 
maximum velocity resolution (1.3 km s^^) showing the Hi emission distribu- 
tion at different velocities were constructed. The Hi emission is shown in Fig. 
2 for HD 38666, Fig. 4 for HD 124979, Fig. 6 for HD 163758, and Fig. 8 for 
HD 171589. Each individual image corresponds to a velocity interval of 4 km 
s~^. Although the whole velocity range at which Hi emission is detected was 
analyzed, the figures include only the images where Hi structures probably 
linked to the stars are detected. 

The identification of an Hi structure associated with a certain star requires 
a careful inspection of the Hi images, looking for shells and voids that might 
be related to the star. Several conditions are necessary to associate an Hi 
structure with a certain star: (a) the star is expected to be located inside 
the Hi void or close to its inner borders; (b) the kinematical distance to 
the structure should coincide, whithin errors, with the spectrophotometric 
distance; and (c) bearing in mind that the velocity dispersion in the interstellar 
medium is about 6 km s~^, an Hi feature should remain detectable for a larger 
velocity interval to be considered as a physical structure. 

We have found a number of Hi structures towards the target stars that 
fulfill these conditions and, consequently, can be related to the stars. In the 
following sections we describe the Hi structures along with their counterparts 
at other wavelengths. 

The main physical parameters of each structure are summarized in Table 
2, which is described in § 6.1. 



5.1. The ISM around HD 38666 

The Hi emission distribution towards this star within the velocity interval 
from -11.8 to +21.1 km is displayed in Fig. 2. HD 38666 is indicated by a 
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Fig. 2. Hi emission distribution towards HD 38666, integrated in steps of 4 km s^^. 
The velocity range and the greyscale (in K) of cacii image are indicated in its upper 
part. The contour Unas are from the minimum value of the grayscale to 12 K, in 
steps of 1 K. HD 38666 is indicated by the star. 
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Fig. 3. Upper panel. Hi emission distribution within tiic velocity interval 7.2 to 
+17.0 km s~^ showing the Hi structure related to HD 38666. The grayscale is from 
2 to 5 K, and the contour lines arc 2.5, 3, 3.5, 4 and 5 K. HD 38666 is indicated by 
a star. Central panel. Overlay of the IR emission at 100 /^m (grayscale) and the 
same Hi contours of the upper panel. The grayscale corresponds to 1.8 to 2.8 MJy 
ster~^). Bottom panel. Overlay of the H- Alpha Composite image (grayscale) and 
the same Hi contours of the upper panel. The grayscale is from 4.3 to 6.5 R. 
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star. The images show that the 0-type star is projected onto the border of a 
low emission region identified within the velocity interval from -7.7 to +17.0 
km s"""^, centered approximately at (l,b) — (235?5,-27?5). Identification of Hi 
structures is difficult for velocities higher than +17 km s~^ due to the lower 
general Hi emission in the region. 

The cavity and surrounding shell around this star is shown in Fig. 3 
(upper panel). The image was obtained by integrating the Hi emission within 
the velocity interval from -7.7 to +14.0 km s^^. The feature is elongated in 
the direction of the stellar proper motion, which is marked in the figure with 
a black arrow. The central panel displays an overlay of the IRIS image at 100 
/Ltm (grayscale) and the Hi image (contour lines). The Hi shell can partially 
be detected in the far infrared. The presence of a strong IR filament at b « 
-25°is compatible with the existence of a thin Hi wall in this region. An IR 
emission band is also projected onto the cavity. If related to the structure, 
this emission might be linked to the approaching or receding caps. 

The bottom panel shows a superposition of the same Hi contour lines 
and the H-Alpha Composite image. Optical emission is present inside the 
cavity and appears bordered by the neutral envelope towards lower galactic 
longitudes and higher negative galactic latitudes. 

Neither catalogued Hii regions nor supernova remnants were found over 
the regions under study. The image at 408 MHz (not shown here) displays 
a ring-like structure centered at (l,b) = (235?0,-28?0), with inner and outer 
semiaxes of 2?lx0?7 and 4?2xl?9, respectively. The radio continuum feature 
is clearly larger than the Hi cavity, and most of it is projected onto the neutral 
envelope, contrary to what is expected for an interstellar bubble or Hii region. 
Thus, the lack of coincidence between the radio continuum feature and the Hi 
cavity casts doubts on its relation to the neutral shell. 

The systemic velocity of the structure, defined as the velocity at which 
the cavity has its larger dimensio ns and is bette r defin ed, is +8 km s^^. 
According to the analytical fit by Brand fc Blitz ( 1993h . the low systemic 



velocity derived for the Hi structure indicates a kinematical distance of about 
1.0 kpc, in agreement with the stellar distance listed in Table 1. The location 
of the star with cavity, together with the agreement between the kinematical 
and spectrophotometric distances strongly suggest that the Hi structure is 
related to HD 38666. We have adopted a distance d = 0.5+0.1 kpc for this 
structure. 

At this distance, the stellar tangential velocity is quite large, vt w 57+15 
km s~^, and the massive star could have created a bow-shock structure. The 
presence of a bow-shock lik e object related to this star was investigated by 
van Buren &: McCrav ( 1988h . who were not succesful in finding such structure 



in the IRAS images. 

5.2. The ISM around HD 124979 



The analysis of the Hi emission distribution maps towards HD 124979 al- 
lows identification of a cavity and shell probably associated with the star. The 
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Fig. 4. Hi emission distribution towards HD 124979, as in Figure 2. For the seven 
first maps, the contour lines are from 2 to 4 in steps of 0.5 K and for 5 to 10 in steps 
of 1 K, and for the last map, from 8 to 22 in steps of 2 K. 
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Fig. 5. Upper panel. Probable Hi structure related to HD 124979. The image shows 
the Hi column density distribution within the velocity interval -47.9, -34.6 km s~^. 
The grayscale is from 3.5 to 10 K, and the contour lines are from 3.5 K to 12 K 
in steps of 1 K. Central panel. Overlay of the far IR emission distribution at 100 
/xm (grayscale from 21 to 35 MJy ster~^) and the Hi column density distribution 
(contour lines). Bottom panel. Overlay of the H- Alpha Composite image (grayscale) 
from 10 to 14 R and the Hi column density distribution (contour lines). 
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Hi minimum is detected within the velocity interval from -52.0 to -23.2 km 

s~^, decreasing in angular size for velocities v > -31.4 km s~^. At the stellar 
position, the contour lines are distorted in the velocity range -52.0 to -39.7 
km s~^. 

The Hi minimum is better defined from -52.0 to -31.4 km s~^, and thus 
the systemic velocity of the structure can be derived as the central value of the 
mentioned interval, Vgys ~ -42±10 km s~^. Taking into account a velocity 
dispersion of 6 km s~^, the systemic velocity corresponds to a kinematical 
distance of 3.0±1.0 kpc, which is compatible with the spcctrophotometric 
distance. We have adopted d — 3.5±1.0 kpc as the distance to the Hi structure. 

Figure 5 (upper panel) shows the Hi brightness temperature distribution 
towards HD 124979, averaged in the velocity interval from -47.9 to -35.6 km 
s^^, for which the Hi cavity is better defined. In order to delineate the cavity 
we have considered the brightness temperature contour corresponding to 3.5 
K. The centroid of the structure is (l,b) = (317?2,+9?3). Both the cavity 
and the almost complete shell are elongated in the direction of the tangential 
motion of the star (marked with a black arrow). 

The 100 /Ltm IRIS image shows far IR emission bordering the region towards 
the galactic plane, revealing the presence of dust associated with the bright 
portions of the surrounding shell (Fig. 5, central panel). The bright IR areas 
detected inside the cavity, which extend well beyond the cavity, are probably 
unrelated to the Hi structure. The proper motion of the star corresponds to 
a tangential velocity of vt « 200±30 km s~^. A careful inspection of the IRIS 
images does not allow to detect a bow-shock like structure at IR wavelengths. 

Neither catalogued Hll regions nor supernova remnants are linked to the 
structure. The bottom panel of Fig. 5 indicates that no significant opti- 
cal emission is detected towards this region. An inspection of the image at 
408 MHz do not show identifiable radio emission probably related to the Hi 
structure. 

No additional OB stars at a compatible distance are projected onto the Hi 
structure. 

5.3. The ISM around HD 163758 

Figure 6 displays Hi brightness temperature images corresponding to the 
velocity interval from -26.3 to -1.6 km s~^ in the vicinity of HD 163758. An 
inspection of the images shows the presence of a low emission region at 6 < 
-5°. 

The Hi cavity centered at (l,b) = (356?0, -6?5) can be appreciated at 
velocities in the range v from -26.3 to -5.7 km s~^. HD 163758, which is 
indicated by the star symbol, appears projected near the higher density border 
of the cavity, indicated by the piled-up of the contour lines. 

The upper panel of Fig. 7 displays the mean Hi brightness temperature 
in the range (-22.2,-5.7) km s~^, where the cavity and shell centered at (l,b) 
= (356?0, -6?5) are clearly identified in contours and grayscale. The central 
and bottom panels show overlays of the Hi emission distribution and the Ha 
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Fig. 6. Hi emission distribution towards HD 163758, as in Figure 2. For the maps 

corresponding to the interval -26.3 to -5.7 km s~^, the contour hnos are from the 
minimum value of the grayscale to 12 K, in steps of 1 K. For the image showing the 
interval -5.7 to -1.6 km , the contour lines are from the minimum value of the 
grayscale to 40 K, in steps of 2.5 K. HD 163658 and WR109 are indicated by the 
star and the cross, respectively. 
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Fig. 7. Upper panel: Hi emission towards HD 163758 integrated within the velocity 
range -22.2 to -5.7 km s~^. The grayscale is from 6 to 12 K. The contour lines 
arc 4.5 and from 5 to 12 K, in steps of 1 K. Central panel: Overlay of the H- Alpha 
Composite image (grayscale from 32 to 45 R) and the Hi emission (contour lines, 
4.5, 5 and from 6 to 12 K, in steps of 2 K). Bottom panel: Overlay of the IR emission 
at 100 ^m (grayscale from 50 to 80 MJy ster~^) and the Hi emission contour lines 
of central panel. 
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Fig. 8. Hi emission distribution towards HD 171589, as in Figure 2. The grayscale 
is from 5 to 30 K. The contour lines are from 5 to 11 K, in steps of 2 K, and from 
15 to 30 K, in steps of 5 K. 
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Fig. 9. Upper panel: Hi column density distribution towards HD 171589, within the 
velocity range +26.3 to +44.8 km s~^. The grayscale is from 7 to 20 K. Contour 
lines are 7, 8, and 10 K, and from 15 to 30 K, in steps of 5 K. HD 171589 is indicated 
by a five pointed star, and the SNR G17.4-2.3 and SNR G17.8 2.6 are marked by 
crosses. Central panel: Overlay of the H- Alpha Composite image (grayscale from 35 
to 70 R) and the Hi column density distribution shown in the upper panel (contour 
lines). Bottom panel: Overlay of the IR emission distribution at 100 fiui (grayscale 
from 70 to 130 MJy ster"'^) and the Hi column density distribution shown in the 
upper panel (contour lines). 
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TABLE 2 

PHYSICAL PARAMETERS OF THE Hi STRUCTURES AROUND THE O STARS 



HD 38666 HD 124979 HD 163758 HD 171589 



(l,b) center 


(235? 5 -27? 5) 


(317?2,+9?3) 


(356? 0,-6? 5) 


(18? 7,-3? 0) 


Velocity interval vi, V2 (km s^^) 


-7,+17 


-53,-25 


-22,-5 


+26,-1-45 


Systemic velocity v^ys (km s^^) 


+8±1 


-42±3 


-10±1 


+34±1 


Expansion velocity Vexp (km s~^) 


11 


15 


10 


11 


Kinematical distance (kpc) 


1.0 


3.0±1.0 


3.3±1.0 


3.0±0.5 


Adopted distance (kpc) 


0.5±0.1 


3.5±1.0 


3.5±1.0 


3.0±0.6 


Radius of the Hi cavity Rcav 


1?6 


1?25 


0?6 


0?6 


Radius of the Hi structure Rs 


3?0 


1?9 


1?2 


1?5 


Radius of the Hi structure Rs (pc) 


26±5 


116±33 


70±20 


52±9 


Mass in the shell (Mq) 


600±240 


58800±33600 


4300±2450 


6900±2300 


Mass deficiency (Mq) 


80±32 


6000±3400 


1100±630 


350±120 


Swept-up mass Ms (Mq) 


340±140 


32400±18500 


2700±1540 


3600±1200 


ne (cm-^) (f=1.0) 




0.5 




0.9 


M. (Mq) (f=1.0) 




24000 




2900 


< (cm-^) (f=0.3) 




0.9 




1.8 


M; (Mq) (f=0.3) 




13200 




1600 


Kinetic energy (10''* erg) 


0.4 


73-103 


2.7 


4.4-6.3 


Dynamical age (lO'' yr) 


1.2 


4.3 


3.5 


3.8 


Ambient density Uo (cm~^) 


0.2 


0.2-0.3 


0.08 


0.25-0.35 



and 100 /im IRIS emissions, respectively. The Ha image reveals a shell-like 
structure of about 1?5 in diameter centered at (l,b) = (356?0, -6?5), with 
the O-type star projected onto one of its borders. The IRIS image shows 
the lack of IR emission in the region of the Hi cavity. Bearing in mind the 
angular resolution of the Hi data, the correlation of the optical feature with 
the small Hi structure is good, suggesting that the Hi feature is the neutral 
gas counterpart of the optical shell. The fact that the star appears projected 
onto the highest density border of the Hi structure reinforces the association 
between the star and the cavity. 

A search for other massive stars projected onto t he region shows that 
the Wolf-Rayet star WR109 (= V617 Sgr, WN5h-H?. I^n der Huchtl boOll) 
is placed about 2° fa r away from the position of the O star. According to 



van der Huchti (j200lh . WR 109 is located at a photometric distance of 34 kpc. 



The faraway estimated distance is the result of the low optical absorption 
derived for the star, along with its low apparent magnitude. This WR star, 
marked by a cross symbol in Figs. 6 and 7, is also near the border of the 
cavity. Its position well outside the optical shell suggests that the WR star is 
unconnected to the Ha shell. 

The systemic velocity of the Hi structure is -10±1 km s^^. The analytical 
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fit to the circular galactic rotation model predicts that material at velocities 
of about -10 km is placed at near and far kinematical distances of 3.3±1.0 
and 13-14 kpc, respectively. The near kinematical distance agrees with the 
spectrophotometric distance to HD 163758 (see Table 1). Note that kinemati- 
cal distances have large uncertainties in this section of the Galaxy. We adopt 
3.5±1.0 kpc as the distance to the Hi feature. 

According to SIMBAD, no Hii regions or SNR have been identified in this 
area. No additional OB stars w ere found projecte d onto the Hi structure. 
The analysis of the CO images bv lDame et al.l (j200ir) did not show detectable 
molecular emission in this region. 



5.4. The ISM around HD 171589 

Figure 8 displays the Hi emission distribution within the velocity interval 
from +20.1 to +57.2 km s^^ in a large area around the massive star. Two 
cavities are clearly detected for velocities in the range +26 to +45 km s^^, 
centered at (l,b) = (18?7 -3?0) and (17?5 -2?7). HD 171589 appears projected 
onto the cavity at larger galactic longitudes. 

The upper panel of Fig. 9 displays the Hi column density distribution 
within the velocity range from +26.3 to +44.8 km s^^. Both cavities are 
separated in this image. The 0-type star is projected onto the Hi hole at (l,b) 
— (18?7,-3?0), partially surrounded by an Hi shell. 

The central and bottom panels of Fig. 9 show the superposition of the Hi 
column density image (in contours) and the Ha and IRIS images, respectively. 
The central panel reveals a region without optical emission coincident with 
the cavity at (l,b) — (18?7,-3?0). The bottom panel shows that both Hi holes 
are partially outlined by bright IR emission at & > -3? 5. The distribution 
of the Hi, optical, and IR emissions is compatible with the presence of an 
interstellar bubble driven by the massive star. 

The systemic velocity of the structure at (l,b) = (18?7,-3?0) is +36±2 km 
s~^, corresponding to near and far kinematical distances of 3.2+0.5 and 13 
kpc, respectively. The near kinematical distance agrees with the spectropho- 
tometric distance derived for the star. Consequently, we adopt d — 3.0+0.5 
kpc. 

Six sup ernova remn ants were detected in this area at galactic latitudes 
b > -2.8 dCreenl l2004h . Two of them (G17.4-2.3 and G17.8-2.6) appear 
projected o nto the cavity centere d at (l,b) — (17?5,-2?7). Based on the S- 
D relation, ICuseinov et al.l (j2003r ) derived distances of about 6.3 kpc to the 
SNRs. As discussed by many authors (e.g. Green 2005) distances derived 
from the E-D relation are hardly reliable. Both SNRs are detected at 35 
MHz ( Dwarakanath fc Udava Shankarill990r ). An inspection of the image at 
4850 MHz shows weak emission also probably related to these remnants. Two 
different facts can be proposed as the origin of the Hi cavity. On one hand, 
the Hi cavity and shell can be associated with the SNRs , as was found in 
many SNRs (e.g. iParon et al.ll2006t iRevnolds et ahlEoOSh . In this case its 
systemic velocity of +36 km s""'^ suggests near and far kinematical distances 
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of 3.2 and 13 kpc, respectively. On the contrary, if the Hi hole originates in 
absorption due to the radio continuum sources, the distance to the SNRs can 
be inferred from the highest positive velocity at which the hole is detected 
(about +52 km s^^). This suggests that the SNRs are placed at distances > 
4.5 kpc. Further data are necessary to elucidate this question. 



CO emission (jPame et al.l 120011 ) is present mainly at & > -3°. The inte- 
grated CO emission (not shown here) in a velocity range similar to that of 
the Hi emission does not show molecular material linked to the hole at (l,b) 
= (17?5,-2?7). Weak CO emission present in the interval +28.0 to +43.6 km 
s~^ encircles the border of the cavity related to the SNRs towards b > -3°. 

6. DISCUSSION 

6.1. Physical parameters 

The main physical parameters of the neutral gas structures linked to the 
0-type stars are summarized in Table 2. The (l,b) centers correspond to 
the approximate centroid of the features. The velocity interval indicates the 
velocity range where the Hi structures can be identified, being vi and V2 the 
lowest and hi ghest velociti e s, resp ectively, at which the features are detected. 



Following ICappa et al.l (|2008l ) , the expansion velocities were estimated as 
Vcxp — (v2 — vi)/2 + 1.3 km s~^. The derived values are lower limits since 
the caps of the expanding shells are not detected in the present cases. The 
extra 1.3 km takes into account the fact that the caps may be present 
just outside the velocity range at which the Hi cavity is detected, that is at 
velocities vi - 1.3 km s^^ and V2 +1.3 km s~^. 

For the four selected stars the adopted distances are compatible with 
the spectrophotometric and kinematical distance estimates. Uncertainties in 
the adopted kinematical distances for HD 124979, HD 163758, and HD 171589 
arise in a velocity dispersion of ±6 km adopted for our Galaxy. 

The radius of each cavity corresponds to the geometric mean of the major 
and minor semiaxes, while the radii of the Hi shells (Rs) were evaluated from 
the position of the maxima in the envelopes. Errors in radii conic from the 
distance uncertainty. 

Neutral atomic masses in the shells and mass deficiencies in the cavities are 
also hsted in Table 2. The structures surrounding HD 124979 and HD 171589 
are easier to define over the halves further away from the galactic plane. In 
each case, the other half is contaminated with diffuse emission from the galac- 
tic plane. Thus, for these two structures, the neutral mass in the shells were 
derived as twice the mass associated with the better defined half. The sewpt- 
up neutral masses were obtained as mean values between the neutral mass 
deficiency in the voids and the neutral mass in the envelopes. This procedure 
allows to remove a first order contribution of the background emission, since 
both mass determinations are probably contaminated with neutral gas unre- 
lated to the structures. Values listed in Table 2 include a typical interstellar 
He abundance of 10%. 
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TABLE 3 

ENERGETICS OF THE STRUCTURES 





HD 38666 


HD 124979 


HD 163758 


HD 171589 




34600 


37100 


40200 


40000 


log L (Lq) 


5.0 


5.6 


6.03 


5.7 


M (Mq) 


23 


40 


70 


47 


M {IQ-^ Mq yr"^) 


0.00032''-0.005'^ 


1.8 ° 


8.3" 


2.7° 




0.02'*-0.35'= 








Voo (km s-i) 


1200'' 


2100* 


2400' 


2500* 


L™ (lO^s erg s-i) 


0.0002-0.16 


2.5 


15.0 


5.4 


Dynamical age (lO'' yr) 


1.2 


4.3 


3.5 


3.8 


E™ (10^« erg) 


0.008-6.0 


350 


1660 


660 


e 


50-0.07 


0.2-0.3 


0.002 


0.007-0.01 


Stellar lifetime*^ (10® yr) 


7.0 


4.5 


3.0 


4.0 



' jMartins etall l|2005l ) 



IChlebowski fc Garmanvl (|l99ll ) 
' iHowarth fc Prinial l|l989ll 
"Estimated from lVink et al.l l|2000l ') 
'' Adopted fromlPrinia et all (|l990l l 
^Schaller et Zl l|l992[ ) 



In the previous sections we have stated that no radio continuum emission 
associated with the Hi structures was detected. Assuming that the massive 
stars have ionized the surrounding gas through their strong UV photon flux, 
and that this ionized material has been swept-up and is present in the inner 
borders of the neutral envelopes, we can derive an upper limit for the electron 
density Ug and the ionized mass from the rms flux density at 4.85 GHz 
(7.7 mJy beam~^). Adopting Rou = 1-0 Rcav and i?i„ — 0.9 Rcav as the 
outer and inner radii of the ionized regions, we estimated upper limits for 
the flux densities S4,s5ghz = 3.7 and 1.0 Jy for the ionized regions around 
HD 124979 and HD 171589, respectively. Estimates of the physical parameters 
of the Hli regions tow a rds th ese stars can be obtained using the expressions by 



Mezger fc HendersonI (jl967l ) for spherical ionized regions of constant density. 
Adopting an electron temperature of 8x10'^ K and a volume filling factor / 
— 1.0, we derived the electron density and the associated ionized mass Mi. 
A different electron density and ionized mass can be estimated by considering 
an alternative filling factor. For an ionized shell with outer and inner radii of 
i?oM and Rim in which the plasma covers an area A equal to 50% of the surface 
of the shell, the filhng factor can be derived as / = Ax {Rl^ - R^n)/ R%u = 
0.3. This value was used to derive and M/. The derived values for are 
consistent with the high z-distances of the stars. Note that electron densities 
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and ionized masses are upper limits. Unfortunately, the lack of radio data at 
frequencies higher than 1 GHz precludes from deriving upper limits for the 
regions of HD 38666 and HD 163758. 

The kinetic energy Ey^ = A/bVcxp^/2 was estimated from the expansion 
velocities listed in Table 2 and the neutral and ionized masses in the structures. 
The range in kinetic energies corresponds to the fact that we have taken into 
account the neutral atomic mass only, and the neutral atomic and ionized 
masses. 

Dynamical ages were estimated as td — 0.55x10^ Kg/vexp yr (|McCrav 



1983I ). where R^ is the radius of the bubble (pc), Voxp is the expansion velocity 



(km s~^), and the constant represents a mean value between the energy and 
the momentum conserving cases. 

Finally, ambient densities Uo were derived by uniformly distributing the 
associated mass within a sphere of radius R^. For the regions of HD 124979 
and HD 171589 two values are listed in Table 2. The first one was obtained by 
distributing the swept-up neutral mass (Ms) while the second one was derived 
by distributing the neutral and ionized masses {Ms + M-). The major source 
of error in radii and masses is the distance uncertainty. The low ambient 
densities are consistent with the structures being far from the galactic plane. 



6.2. Energetics and origin of the structures 

We will analyze here whether the massive stars can provide the energy 
to create the cavities and shells found in the previous sections through their 
stellar winds and UV photon fluxes. To test the former possibility, we will 
estimate the mechanical energy E^, released into the ISM for the massive stars 
and compare it with the kinetic energy i?k of the structures. 

Table 3 summarizes relevant parameters useful to evaluate the energetics of 
the structures. The values of effective temperature T^ff, stellar luminosity L, 
and stellar mass listed in the first three rows correspond to the spectral classi- 
fication and luminosity class of the 0-type stars listed in Table 1. These data 



were u sed to derive the mass loss rate M following the recipe bv lVink et al 



([20001). Based on the mass loss rates and terminal velocities, we estimated 
the mechanical wind luminosity for each star as Lw = 

The stellar wind mechanical energy E^{— L^t) released by the massive 
stars were derived from the dynamical ages of the bubbles. The resulting 
values are included in Table 3. 

The ratio e between the kinetic energy i?k and the mechanical energy E^ 
provided by the massive stars during the dynamical age of the structures is also 
included in Table 3. Evolutionary models of stellar wind bubbles predict that e 
0.2 for the energy con serving case and e < 0.1 for the momentum conserving 



case (see McCrav 19831 ). Although the uncertainty in this value is large (at 



least 70% adopting a 30% error in the distance), it is clear that HD 163758 
and HD 171589 are capable of blowing the observed estructures. The result 
for HD 124979 is still consistent with an interstellar bubble interpretation. 
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The obtained dynamical ages for HD 124979, HD 163758, and HD 171589 
are compatible with t he ages derived from evolutionary tracks for stars with 



solar abundances (see Schaller et al. 1992h . 



The case for HD 38666 is more complex. The large uncertainty in the es- 
timated mechanical energy precludes for giving a clear conclusion about the 
origin of this structure . The derived dyna mical age i s lower than the age 



However, 



Martins et al 



(200. 



i 



estima te an age < (2-4) x 10^ y r, closer to the dynamical age. On the other 



hand, iHoogerwerf et al.l (|200lf ) proposed a binary-binary collision between 
three stars of Trapezium cluster: AE Aur, HD 38666, and l Ori , becoming 
HD 38666 and AE Aur runaway stars. The dynamical ejection scenario took 
place 2.5 X 10^ yr ago. We note that a search for other catalogued massive stars 
at a distance compatible with that of the structure gave negative results. A 
different analysis can be done bearing in mind the observed proper motions. 
HD 38666 and HD 124979 have large spatial velocities, about 57±11 km s^^ 
for HD 38666, and 200±30 km s'^ for HD 124979. Considering these veloci- 
ties, which include the uncertainty in the adopted distance, it took HD 38666 
about (0.4-0.9) X 10^ yr to cross the Hi structure. As regards HD 124979, the 
time necessary to cross the structure was about (0.6-1.0) x 10^ yr. Adopt- 
ing 0.7x10*^ yr for HD 38666 and 0.8x10^ yr for HD 124979, the mechanical 
energy turns out to be (0.004-3.5) x lO-^^ erg for HD 38666 and 65x10"^® erg 
for HD 124979. The large spatial velocity of the runaway star HD 124979 
shortenes the crossing time through the structure, which makes their origin 
to remain unclear. The resulting e-values indicate that other energy sources 
are necessary to create both structures through the stellar wind mechanism. 
Additional studies are necessary to identify the agents. 

No bow-shock like features were found related to HD 38666 and HD 124979. 



As shown bv lRaga et al.l (|l997l ) and lHuthoff fc Kapeij (j200ir ). bow-shocks ap- 
pear associated with only about 30-40% of OB runaway stars. As indicated 
by the last authors, large separation from the galactic plane, extremely high 
space velocities and large distances make difficult the formation and detec- 
tion of bow-shock structures. Moreover, the physical conditions of the ambi- 
ent medium where the star is inmersed play a major role in determining the 
existence of bow-shocks. Particularly, these authors find an anticorrelation 
between bow-shocks and hot bubbles. Whether the observed structures orig- 
inate in the mass flow of the massive stars, the lack of bow-shocks in these 
cases is consistent with their conclusions. 



Taking into account the spectral classification and the luminosity class of 
the associated massive stars and the ambient density where the structures are 
evolving, the estimated radius of the Stromgren's spheres is in all cases larger 
than the Hi cavities. The smaller size of the cavities can be explained taking 
into account that a certain amount of UV photons are used in dust heating 
or drain from the patchy envelopes. 
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7. SUMMARY 

We have analyzed the interstellar medium in the environs of the 0-type 
stars HD 38666, HD 124979, HD 163758, and HD 171589. The study of the 
neutral hydrogen distribution in direction to these stars allowed us to disclose 
Hi structures located at kinematical distances compatible with the stellar 
distances and probably related to the stars. 

Assuming a stellar wind mechanism, the derived dynamical ages for the Hi 
structures related to HD 124979, HD 163758, and HD 171589 are compatible 
with the lifetimes of the 0-type stars in the main sequence, reinforcing the 
association with the stars. 

We have investigated the counterparts of the structures at other frequen- 
cies. The surrounding Hi shells around HD 38666, HD 124979, and HD 171589 
have IR counterparts, revealing the presence of dust associated with the Hi 
gas. The Hi shell related to HD 163758 has IR and Ha counterparts, with the 
star projected onto one of its higher density borders. 

To investigate the origin of the structures, we have compared the mechani- 
cal energy released into the ISM for the massive stars and the kinetic energy of 
the structures. Our results for the energy conversion efficiency for HD 163758 
and HD 171589 indicate that the stars are capable of blowing the observed 
structures. As regards HD 38666 and HD 124979, additional energy sources 
are probably necesary, taking into account the derived e-values and the large 
tangential velocities. Aditional studies are needed to clarify the origin. 
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